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ABSTRACT 

We study the structural characteristic of the variable DA white dwarf G117B-15A 
by applying the methods of asteroseismology. For such a purpose, we construct white 
dwarf evolutionary models considering a detailed and up - to - date physical description 
as well as several processes responsible for the occurrence of element diffusion. We 
have considered several thickness for the outermost hydrogen layer, whereas for the 
inner helium-, carbon- and oxygen- rich layers we considered realistic profiles predicted 
by calculations of the white dwarf progenitor evolution. The stellar masses we have 
analysed cover the mass range of 0.50 < M^/Af© < 0.60. 

The evolution of each of the considered model sequences were followed down to 
very low effective temperature; in particular, from 12500K on we computed the dipo- 
lar, linear, adiabatic oscillations with radial order k = 1, • ■ • ,4. We find that aster- 
oseismological results are not univocal regarding mode identification for the case of 
G117B-15A. However, our asteroseismological results are compatible with spectroscop- 
ical data only if the observed periods of 215.2, 271.0 and 304.4 s are due to dipolar 
modes with k = 2,3,4 respectively. 

Our calculations indicate that the best fit to the observed period pattern of G117B- 
15A corresponds to a DA white dwarf structure with a stellar mass of 0.525 Mq, with a 
hydrogen mass fraction log (Mh/M*) ^-3.83 at an effective temperature T^s f»11800K. 
The value of the stellar mass is consistent with that obtained spectroscopically by 
Koester & Allard. 

Key words: stars: evolution - stars: interiors - stars: white dwarfs - stars: oscillations- 
stars: individual (G117B-15A) 



1 INTRODUCTION 

Asteroseismology has become a powerful method to disen- 
tangle the internal structure and evolution of stars by means 
of the study of their oscillatory pattern. This technique, very 
sophisticated in the case of the Sun, has also undergone a 
strong development in other stars, in particular in variable 
white dwarf (WD) stars (see e.g. Brown & Gilliland 1994; 
Gautschy & Saio 1995; 1996). 
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Pulsating WDs show multiperiodic luminosity varia- 
tions in three ranges of effective temperatures (Tefj ) corre- 
sponding to the currently called DOV, DBV and DAV (see, 
e.g., the review by Winget 1988). Of interest in this work 
are the DAVs (hydrogen - dominated atmospheres), or ZZ 
Ceti, that pulsate in the instability strip corresponding to 
11000 K < Tcff < 13000 K. The periodicities in the light 
curves of pulsating WDs are naturally explained in terms of 
nonradial g- modes of low harmonic degree {£ < 2), driven 
by the "k mechanism" working in a partial ionization region 
near to the stellar surface (Dolez & Vauclair 1981, Winget 
et al. 1982). The periods (P) are found within a range of 
100s < P < 1200s and photometric amplitudes reach up to 
0.30 magnitudes. 

Most asteroseismological studies performed on WDs to 
date rely on stellar models constructed under some simpli- 
fying hypotheses. One of the most relevant is related to 
the profile of the internal chemical composition in the in- 
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terface zones. In such regions, the equilibrium diffusion in 
the trace element approximation has been widely employed 
to infer the profile of the chemical distribution (see Tas- 
soul, Fontaine & Winget 1990). The main motivation for 
considering this approximation is to avoid the solution of 
time dependent element diffusion as the WD evolves. In the 
frame of such an approximation, the profile of the inter- 
face region is very simple: its functional form is a power 
law. The transition zone is separated into two parts: an up- 
per one in which one element is dominant and the other is 
considered as a trace and a lower region in which the role 
of the respective elements is reversed. Because these two 
power law solutions are matched for fulfilling the condition 
of conservation of mass of each element, a discontinuity in 
the derivative occurs just at the matching point. The ex- 
ponent of the power law solution is directly related to the 
state of ionization of the stellar plasma. Thus, the structure 
of the whole interface zone, in the frame of this standard 
treatment, can be modified only if the plasma suffers from 
a modification in the state of ionization as a result of stellar 
evolution. Calculations of asteroseismology of WDs in the 
frame of such a standard treatment for the chemical inter- 
faces are those of, e.g., Bradley (1996, 1998, 2001), Bradley 
& Winget (1994), Brassard et al. (1991, 1992ab), Fontaine et 
al. (1992), Montgomery & Winget (1999), Metcalfe, Nather 
& Winget (2000), Metcalfe, Winget & Charbonneau (2001), 
Montgomery, Metcalfe & Winget (2001). 

From an evolutionary point of view, the shape of the 
chemical interfaces may not be criticaj^ but in pulsational 
studies they provide a non - negligible contribution to the 
shape of the Ledoux term of the Brunt- Vaisala frequency 
(Brassard et al. 1991). Thus, we can, in principle, expect 
differences between studies with equilibrium diffusion and 
those in which other more physically sound treatment is 
performed. This is a very important point to be made in 
connection with the aim of the present work. 

Since the pioneering work of Iben & MacDonald (1985), 
we know that element diffusion modifies the chemical abun- 
dance distribution within a WD star even during evolution- 
ary stages corresponding to the ZZ Ceti domain (see Iben 
& MacDonald 1985, particularly their Fig. 4). Few calcula- 
tions exist in the literature in which the evolution of WDs 
is addressed in a self consistent way with time dependent 
element diffusion. Amongst them, we mention the study 
of MacDonald, Hernanz & Jose (1998) aimed at studying 
the carbon pollution in cool WDs. Also, Dehner & Kawaler 
(1995) used non - equilibrium diffusion together with evo- 
lutionary calculations to study WDs with helium envelopes. 
Finally, Althaus, Serenelli & Benvenuto (2001) have recently 
shown that diffusion induces the occurrence of thermonu- 
clear fiashes in helium - core WDs, causing the evolution 
of such kind of WDs to occur on timescales significatively 
shorter than predicted by models without diffusion. This 
has been particularly important in solving the discrepancy 
(Van Kerkwijk et al. 2000) about the age of binary systems 
containing a millisecond pulsar and a helium WD. 

As far as we are aware, the only work aimed at ex- 
ploring the role played by diffusion in the period and its 

^ An obvious and important exception is when the tail of chem- 
ical profiles are subject to nuclear burning. 



rate of change for g-modes in pulsating DA WDs is that of 
Corsico et al. (2001b). In particular, the authors found that 
the differences in the shape of the chemical profiles at the 
interface zones induce appreciable changes in the periods, as 
compared to the case of equilibrium diffusion in the trace el- 
ement approximation. Also, there are noticeable changes in 
the period derivative, which are due in part to the evolution 
of the chemical profile during the cooling of the WD across 
the ZZ Ceti instability strip. 

Since sometime ago there have been several works avail- 
able in the literature in which the observed period structure 
of a particular object is fitted to theoretical predictions. Such 
works show that, in principle, information about the stellar 
mass and the stratified outer layer structure can be inferred 
(see, e.g., Winget et al. 1991; Fontaine et al. 1992; Pfeiffer 
et al. 1996; Bradley & Kleinman 1997; Bradley 1998, 2001; 
Bradley & Winget 1994). As these studies are based on a 
simplified treatment of diffusion, we believe that it is worth 
revisiting this problem on the basis of more detailed models 
of WD structure. In particular, in view of the results found 
in Corsico et al. (2001b), we expect to find differences in 
the period fitting to a particular object, as compared with 
the situation in which the standard treatment is used. The 
only way to find how important can be such differences is 
by performing a detailed asteroseismological study choosing 
a well studied object. 

We consider G117-B15A as an optimal target for our 
study. G117-B15A is an otherwise typical DA WD, the vari- 
ability of which was discovered by Mc Craw & Robinson 
(1976) and, since then, it has been monitored continuously. 
The mass and in particular the effective temperature of this 
star have been the subject of numerous spectroscopic rede- 
terminations (see Gautschy, Ludwig & Freytag 1996 for a 
summary), in particular values of 0.59 M0and 11620 K, re- 
spectively, have been derived by Bergeron et al. (1995). More 
recently, Koester & AUard (2000) (hereafter KA) have sug- 
gested a somewhat lower value for the mass of 0.53 M0and 
a higher effective temperature of T^s =11900 ± 140 K. 

G117-B15A has periods of oscillation of 215.2, 271 and 
304.4 s (Kepler et al. 1982). Notably, for the 215.2 s mode it 
has been possible to find a value of its temporal derivative 
(Kepler et al. 2000) to be P = (2.3 ± 1.4) x lO"^'^ s s'^ 
Interestingly, the 215.2 s mode present in G117-B15A is the 
most stable optically registered oscillation with a stability 
comparable to that of the most stable millisecond pulsars. 
As G117-B15A is a well know oscillator, it has motivated 
the interest of several researchers in its internal structure. 
For example, Bradley (1998) has found the best fitting to 
the period pattern with a model of ~ 0.6 MQand depending 
on the identification of the modes (see below for an expla- 
nation of its meaning), the favoured value for the hydrogen 
envelope mass fraction Mh/M, is -7 < log(MH/M,) < —6 
or -5 <, log (Mh/M.) <, -4, and log(MHc/M*) =-2. More 
recently, Corsico et al. (2001a) have found, on the basis of 
WD models in which diffusion is neglected, the best fitting 
to the period pattern with a 0.55 MqWD model with a car- 
bon - oxygen interior, log(MHc/M,) =-2 and log(MH/M,) 
=-4 as predicted by stellar evolution. Notably, such a fitting 
is in nice agreement with one of those proposed by Bradley 
(1998). 

It is the aim of this work to perform a fitting to the pe- 
riod structure present in G117-B15A by computing the non- 
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radial eigenmodes in the frame of linear and adiabatic ap- 
proximation, and evolutionary WD models in which time de- 
pendent element diffusion is properly accounted for. About 
the way we shall perform such a fitting, some words are in 
order. In handling models like those we shall employ it is 
quite obvious that we are not in a position to employ tech- 
niques like those used in the so - called genetic algorithm 
(see its application for the case of the DBV GD358 in Met- 
calfe et al. 2000). In the context of ZZ Ceti stars a similar 
approach though to a less extent has been repeatedly em- 
ployed by numerous investigations (e.g the complete and de- 
tailed study performed by Bradley 1998, 2001). Up to now, 
the shape and thickness of the relevant chemical interfaces 
are usually treated as a free parameter. However, a more 
physically sound treatment can be performed when account 
is made of time dependent element diffusion in evolutionary 
models of WDs. This aspect of WD evolution is one of the 
most important when an attempt is made to compare obser- 
vations with theoretical expectations from pulsating WDs. 
Unfortunately, there are currently some uncertainties (e.g., 
the treatment of convection, the rate of the critical nuclear 
reaction ^^C(q:, 7)^®0, wind mass loss episodes, etc.) that 
prevent us, provided a mass value and pre- WD composition, 
from predicting a whole definite internal structure. This is 
especially true for Mh /M« . Thus, in view of these facts, we 
shall consider the value of Mh/M, as a free parameter. We 
assumed log (Mhc/M*) —-2 throughout this paper. 

The remainder of the paper is organized as follows: In 
Section ^ we discuss the physical ingredients we have em- 
ployed together with the computational strategy we have 
employed. Section ^ is devoted to presenting the details of 
the evolutionary models we have constructed down to the 
conditions relevant for the WD we are investigating. In sec- 
tion ^ we present the asteroseismological analysis we per- 
formed for the particular case of G117B-15A. Section ^ is 
devoted to the discussion of the results we have found and 
finally in Section ^ we give some concluding remarks and 
insights for future investigations. 



2 DETAILS OF THE COMPUTATIONS 

For this work we have employed the same code and physical 
ingredients as in Corsico & Benvenuto (2002) and Corsico 
et al. (2001ab). In particular for the diffusion processes we 
considered gravitational settling, chemical and thermal dif- 
fusion as in Althaus & Benvenuto (2000). 

In our code, evolutionary and pulsational calculations 
are performed in a fairly automatic way. After selecting a 
starting stellar model we choose an interval in P and Tcff . 
The evolutionary code computes the model cooling until the 
hot edge of the T^s -interval is reached. Then, the program 
calls the set of pulsation routines beginning the scan for 
modes. When a mode is found, the code generates an ap- 
proximate solution which is iteratively improved to conver- 
gence (of the eigenvalue and the eigenmode simultaneously) 
and stored. This procedure is repeated until the period in- 
terval is covered. Then, the evolutionary code generates the 
next stellar model and calls pulsation routines again. Now, 
the previously stored modes are taken as initial approxima- 
tion to the modes of the subsequent model and iterated to 
convergence. Such a procedure is automatically repeated for 



all evolutionary models inside the chosen Tcff - interval. The 
computational strategy described above has been success- 
fully applied in fitting the observed periods of G117-B15A 
to impose constraints on the mass of axions (Corsico et al. 
2001a) in the frame of WD models less detailed than those 
we shall employ in the present work. We refer the reader 
to Corsico & Benvenuto (2002) and Corsico et al. (2001ab) 
and references therein for further details. Let us quote that 
most of our evolutionary models have been divided in more 
than 2000 meshpoints, whereas for mode calculations we em- 
ployed 5000 meshpoints. 

In order to start our evolutionary calculation we have 
employed an artificial heating technique detailed in Althaus 
& Benvenuto (2000). Let us quote that in this technique, the 
stellar model is forced to undergo some unphysical evolution 
by introducing an artificial energy release that, when the 
star gets bright enough, is switched off smoothly. Then, the 
star relaxes to the physical cooling branch after few tens of 
models. 

In employing this artificial technique, we must be very 
careful in avoiding generating an artificial chemical profile. 
Thus, in order to verify that the star not only relaxes to the 
physical cooling branch but simultaneously relaxes to a cor- 
rect internal chemical profile, we have employed, for a fixed 
log (Mh/M«) value, different shapes of the initial hydrogen 
and helium profiles. We found that, if the initial model is 
far away enough from the evolutionary stage we are inter- 
ested in, the fine details of the initial profile are irrelevant 
because diffusion largely evolves it to a well defined one. As 
a matter of fact, the structure and the chemical profiles are 
physically plausible far before the star attains the effective 
temperature range relevant for G117-B15A. 

A better approach for constructing pulsating WD mod- 
els would obviously be to start computing evolutionary mod- 
els from much earlier in the star history. This possibility, 
better from a physical point of view, again would require 
an enormous amount of computing time for a definite se- 
quence. Thus trying to fit a period pattern in such a way 
is, in our opinion, far beyond what seems reasonable at the 
present state of the art. In any case, it would be interest- 
ing to explore the pulsational properties of models resulting 
from a full evolutionary calculation. We plan to present such 
a study in a future paper. 

To be specific, we have started the calculations with 
models of log(L/L0) ~ 2.05, log Teg fti4.46 and masses of 
0.50, 0.55, and 0.60 AIq (in performing the fitting of the mod- 
els to observations we have also computed supplementary 
models of 0.525 and 0.5375 Mq). The mass fraction em- 
braced by the outermost hydrogen layer has been assumed in 
a way that, after evolving along the cooling branch down to 
Tcff =12500 K the values of log (Mh/M.) present in the star 
are those given in Table 1. At the starting conditions we have 
assumed the Salaris et al. (1997) core composition, whilst 
for the outer layers the CNO abundances were those corre- 
sponding to the pre- WD evolution. In order to get different 
values of Mh/M, , we have started from the model with the 
largest Mh/M» and simply replaced by ''He at the base 
of the hydrogen envelope. Let us quote that the initial model 
has internal temperatures low enough that prevent the oc- 
currence of nuclear reactions up to the moment at which the 
WD model has relaxed to the physical branch. We should 
mention that models become physically adequate (after the 
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unphysical transitory) for describing a WD at log (L/L©) ~ 
1.4, logTeff ~4.83, far from the conditions relevant for DAV 
WDs. Also, the values of log(MH/M,) at the instability 
strip are somewhat lower than those included in Table 1 
simply because nuclear reactions are still operating at such 
evolutionary stages. However this effect is noticeable only 
for the case of the thickest hydrogen envelopes on the most 
massive objects considered here. 

Table 1. Values of log{MH/M*) for our computed models 
at Tcff =12,500 K 

O.SOA/q O.55M0 O.6OM0 O.525M0 O.5375M0 

-3.815 -3.862 -3.941 -3.831 -3.843 

-4.193 -4.224 -4.249 -4.198 

-4.679 -4.684 -4.692 

-5.175 -5.177 -5.180 

-5.671 -5.671 -5.672 

-6.160 -6.158 -6.156 -6.159 

-6.640 -6.633 -6.700 -6.634 

-7.071 -7.047 -7.028 

-7.423 -7.385 -7.349 



To our notice, this is the first opportunity in which car- 
bon - oxygen, intermediate mass WD models with a set of 
values in Mh/M, are evolved in the frame of time dependent 
element diffusion. Prior to this work, the kind of calcula- 
tions carried out here were performed only for the values of 
Mh/M, predicted by stellar evolution theory (Iben & Mac- 
Donald 1985; 1986). 



3 EVOLUTIONARY RESULTS 

For each initial model, we have computed its evolution down 
to the conditions corresponding to the DAV instability strip. 
In doing so, we computed about ten thousand models for 
each sequence. Perhaps the most relevant characteristic of 
the models for the kind of study we are performing here 
is the internal chemical profile. Results representative for 
the case of the models corresponding to the DAV instabil- 
ity strip are presented in Fig. 1. In this figure we depict 
the internal chemical profiles for a 0.55 MqWD model at 
log Toft ^ 4.096 for log(MH/M.) =-3.862. The adopted ef- 
fective temperature is slightly higher than the last estima- 
tion for G117B-15A (logTeff « 4.076). Thus, these profiles 
should be representative of the conditions expected for the 
interior of G117B-15A in the frame of our set of models. 

Let us discuss the main characteristics of these profiles. 
Before discussing them, we have to warn the reader that the 
way we have modified the value of Mh/M, is ad-hoc. In fact, 
it is not clear at present if there exists, for a fixed stellar mass 
value, an evolutionary mechanism capable to produce a dis- 
tribution of Mh/M« values differing each other by orders of 
magnitude. In addition, the outward tail of the distribution 
of initial heavy element abundances may be different from 
the results we shall present below. However, we should re- 
mark that the profile of the light elements present in the 
outer parts of the models should be essentially correct. 

With regard to the chemical profiles of our models we 
want to mention that for log(l — Mr/M,) > —1.6 the car- 
bon - oxygen distribution is that predicted by pre- WD evo- 



lutionary calculations and it is the same for all of our evolu- 
tionary sequences. However this is not the case for the chem- 
ical structure of the outer layers, which are, in turn, crucial 
for the pulsational properties of the models. The main vari- 
ations are related to the whole distribution of ^H, ^He and 
the outwards tail of *He. While a variation in the assumed 
value of Mh/M, obviously implies a modification in the pro- 
files, we find a very noticeable variation in the slope of the 
profiles for the isotopes present in the outer layers of our 
WD models. In the model with log (Mh/M,) =-3.815 it is 
remarkable the large slope in the tail of the hydrogen dis- 
tribution. In spite of this steep profile, degeneracy prevents 
hydrogen from diffusing further inwards. In models with a 
smaller hydrogen content, the tail of the hydrogen distri- 
bution is smoother. Note also that the shape of the tail of 
the hydrogen distribution at the evolutionary conditions rel- 
evant for DAV WDs is strongly dependent on the value of 
the assumed Mh/M» . 

Regarding the tail of the ''He distribution, this is largely 
due to gravitational settling that forces ''He to sink down. 
In particular our models also present a zone in which the 
chemical composition is largely dominated by helium. Also, 
it is of interest the distribution of ^He. It is clearly noticeable 
that most of the ^He is located in stellar layers at which 'H 
and ''He are present. This is a consequence of the interplay 
between nuclear reactions that form "^He by means of the 
first two reactions of the PPI chain (Clayton 1968) 

^H + 'H^^'He 

and gravitational settling. Notice that the outer tail of 
the "^He distribution is similar to the corresponding to ''He 
which is, in part, due to the operation of the reaction 
^He ^He ^ ''He 2 'H. 

As the chemical profile of the hydrogen in surface lay- 
ers is one of the key ingredients in determining the g-modes 
of DAV WDs, we show in Fig. 2 the base of the hydrogen 
envelopes for our suite of 0.55 MQmodels at the ZZ Ceti do- 
main. In addition, we show the profiles predicted by the stan- 
dard treatment of equilibrium diffusion in the trace element 
approximation for the same stellar mass and Mh/M« values. 
It is clear that time dependent element diffusion predicts hy- 
drogen profiles that are somewhat different from those given 
by the standard treatment. In the standard treatment, the 
slope of the profile is determined solely by the state of ion- 
ization of the plasma (see Eqs. 26 - 29 of Tassoul et al. 
1990). Thus, if the state of ionization is the same, the slope 
should remain essentially unaltered with further evolution. 
We have verified that, even in the case of the thinnest hydro- 
gen envelopes we have considered, the base of the hydrogen 
envelope is located at layers with temperatures above lO^'K, 
which ensures complete ionization for hydrogen and helium 
and this is so even at the effective temperatures correspond- 
ing to the DAV instability strip. Also, notice that the sudden 
change of the slope at the matching point of the two solu- 
tions of the standard treatment (denoted by circles along 
the curves) has a direct consequence on the enhancement of 
the trapping of some modes in the hydrogen envelope (see, 
e.g.. Brassard et al. 1992a, Bradley 1996). Thus, our results 
show the necessity of allowing for a more detailed treatment 
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of the chemical evolution at performing asteroseismology of 
WDs. 

Now, let us discuss the impact of our chemical profiles 
on the Brunt- Vaisala frequency, a key quantity entering as 
a coefficient in the equations of non-radial pulsations. This 
is defined as (Unno et al. 1989) 



,,2 / 1 dlnP d\np 

N = g - 

^ VTi dr dr 



(1) 



where all the symbols have their usual meaning^ and Fi 
is the first adiabatic exponent. However, in order to avoid 
numerical noise, we use the expression given by Brassard et 
al. (1991): 



^2 ^ 9_£XT^y^^^y^Q^ 

P Xp 



(2) 



The Ledoux term _B, for the case of a M-component 
plasma, is given by 



B 



M-l 

1 ■^-^ d In Xi 
XT ^ dlnP 

i—l 



(3) 



where 



Xp = 



XT = 



dlnP 
din p 



/ d\nP \ 
[dlnTj 



T,{Xi} 



P,iX,} 



^ / ainP X 
\dlnXj 



p,T.{Xj^,} 



In Fig. 3 we show the squared Brunt- Vaisala fre- 
quency and the Ledoux term for models of 0.55 A^Qwith 
logTeff «4.096 and all the considered values of Mh/M« . 
We also show in an inset the Brunt- Vaisala frequency at 
the hydrogen - helium interface for the case of the thick- 
est hydrogen envelope according to our computations and 
to the prediction of equilibrium diffusion (thick and thin 
lines, respectively). In order to show the correspondence of 
the structure of these functions with the internal chemical 
profiles of the star we have also included the curves corre- 
sponding to the abundances of "^He and ^'^C. 

The Ledoux term B has a structure that remains largely 
invariant for the stellar core under changes in Mh/M* . From 
the centre outwards, the three peaks are due to steep slopes 
in the carbon profile. The first two are due to the structure 
of the carbon - oxygen interface whereas the other is due 
to the helium - carbon interface. In contrast to the internal 
behaviour of the Ledoux term, the other peak, due to the 
hydrogen - helium transition, is largely modified, not only in 
position but also in its height. Notice that the lower Mh/M* 
the thinner the peak (see Brassard et al. 1991 and Tassoul 
et al. 1990 for comparison of the shape of the B term and 
Bradley 1996 for the profile of the Brunt- Vaisala frequency) . 

It is worthwhile to comment on the fact that because of 



I In this Section P stands for the pressure but in the rest of the 
paper it stands for the period of oscillation. 



our realistic chemical profiles, we would expect to find differ- 
ences in the mode trapping properties as compared with pre- 
vious studies. This important issue surely deserves a careful 
exploration, which will be reported in a subsequent commu- 
nication. 



4 ASTEROSEISMOLOGICAL RESULTS 

An important aspect for the present work is the so-called 
mode identification. That is, the identification of the £ and 
k values corresponding to each observed period. In the case 
of G117-B15A, the results of Robinson et al. (1995) indi- 
cate that the 215.2 s period corresponds to a dipolar mode 
{i — 1). Following the work by Bradley (1998) we shall as- 
sume that the other two modes cited above are also dipolar 
(Brassard et al. 1993; Fontaine & Brassard 1994), and also 
that the other periodicities present in the light curve of the 
star are not associated with actual eigenmodes but are ac- 
tually due to non-linear effects in the envelope (Brassard 
et al. 1993). With regard to the radial order of the modes, 
there exist two possible identifications: following Clemens 
(1994) the observed periods are dipolar modes with k — 2 
(215.2 s), fc = 3 (271 s) and fc = 4 (304.4 s). On the 
other hand, Fontaine, Brassard & Wesemael (1994) iden- 
tified these modes as dipolar with fc = 1 (215.2 s), fc = 2 
(271 s) and fc = 3 (304.4 s). According to these facts, in 
computing the g-modes of our WD models, we have consid- 
ered, as previously mentioned, only £ = 1 and radial order 
fc = 1, • • • , 4. Higher radial orders correspond, for our stellar 
mass values, to periods long in excess compared to those 
observed and thus are not considered here. 

In order to look for a fit to observations we shall consider 
a function ip defined as 



V' = log ^ 



PP 



(4) 



where P^ and Pp axe the i-th observed and calculated pe- 
riods respectively. For the case of G117B-15A, P°= 215.2, 
271 and 304.4 s. As a result of the pulsational calculations 
we have the values of the periods for k = 1, ■ • ■ , 4. Then, 
assuming the two reasonable possible identifications of the 
observed modes (i.e. Pp corresponding to fc = 1,2,3, here- 
after case A, and Pp corresponding to fc = 2, 3, 4, hereafter 
case B) we can compute the function tp in terms of the ef- 
fective temperature of the star for each of the considered 
sequences. Obviously in the case that the fitting were exact, 
ip ^ —oo. Thus the criterium we shall assume in selecting 
the best model for representing G117B-15A is simply to look 
for the minima of the i/) function considering the whole set 
of models. In such a way we shall be able to determine not 
only the internal structure and mass of G117B-15A, but also 
its asteroseismological effective temperature. 

In Fig. 4 we show i/) as a function of the effective tem- 
perature for the case of the set of models with 0.50 Mq. 
In this figure the left panel labeled as A corresponds to the 
identification of the observed modes with fc = 1, 2, 3 whereas 
the right panel B corresponds to the other identification of 
modes we shall consider: fc = 2, 3, 4 modes. As it has been 
found by Bradley (1998), for cases A and B the best fitting 
to observations are found for models with values of Mh/M, 
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very different. While assuming case A we find the best fit 
with rather thin hydrogen envelopes ( log (Mh/M,) = -6.16) 
and in the case B the best fit is found for the case of the 
thickest considered hydrogen envelope: log(MH/M») = - 
3.81. It is remarkable that the minimum found for case A 
is a very deep one, corresponding to a mass and effective 
temperature (Toff «11400K) rather lower than those corre- 
sponding to the last determination of these parameters for 
G117B-15A. 

If we increase the stellar mass to 0.55 Mq , the behaviour 
of the function ip is largely changed (see Fig. 5). Despite that 
in case A the best fitting is for log(MH/M») = -6.63, it 
corresponds however to an effective temperatures far lower 
than that suggested by KA (and outside the DAV instability 
strip). By contrast, in case B the best fit is again provided 
by the model with the thickest hydrogen envelope, but again 
it occurs at effective temperatures somewhat too low. 

Finally, in Fig. 6 we show the behaviour of function tp 
for the case of the models constructed assuming a stellar 
mass value of 0.60 Mq. It is clear that in both cases of 
mode identifications the fits are far poorer than those for 
the other, lower mass values. Thus, the results presented up 
to this point strongly suggest that the asteroseismological 
mass of G117B-15A should be in between 0.50 MQand 0.55 
Mq. Moreover, the only viable mass for the case A seems to 
be 0.50 Mq, because larger mass values give rise to poorer 
fits. 

Let us analyse the results corresponding to case B more 
carefully, centering our attention in models with the thickest 
hydrogen envelopes. From Fig. 4b, it can be noticed that in 
the case of 0.50 i\f0model, the fit is better for increasing ef- 
fective temperatures. However, notice that (Fig. 5b) for the 
0.55 Mq model the situation is just the opposite: the fit im- 
proves for lower effective temperatures. Guided by this fact, 
it can be expected a minimum in tp for some intermediate 
value of the stellar mass. In order to explore this possibil- 
ity we computed models with 0.525 Mq restricting ourselves 
to values of Mh/M* similar to those that provided good 
fits (see Table 1). For this mass value, the identification of 
modes in case A shows a good fit for log (Mh/M.) = -6.63 at 
Tes « 11400 while case B presents an interesting minimum 
at stellar mass and effective temperature values very similar 
to those determined by KA for this object by spectroscopic 
analysis. Finally, in order to perform a finer analysis, we 
compute a sequence of models with 0.5375 MQand a thick 
hydrogen envelope (see Table 1). The results for the whole 
set of masses and thickest hydrogen envelopes for case B is 
shown in Fig. 7. The labels of each curve correspond to the 
values of log (Mh/M.) presented in Table 1. 

The results presented in this Section indicate that, if we 
expect the asteroseismological masses and effective temper- 
atures to be similar to those determined spectroscopically, 
the best fit is found in the frame of case B for a model with 
a stellar mass of 0.525 Mq. For the sake of completeness, we 
also show in Fig. 8 the periods of oscillation for radial order 
A; = 2, 3, 4 for the case of models with the thickest hydrogen 
envelope for each considered value of the stellar mass. 



5 DISCUSSION 

From the results presented in the preceding Section, we con- 
clude that the best fit to observations is found with stel- 
lar mass and Mh/M, values that are very sensitive to the 
adopted mode identification. This has been previously found 
by Bradley (1998). If we restrict ourselves to analyse the 
case of G117B-15A only in the frame of asteroseismological 
models we find no way to be conclusive about the actual val- 
ues for its stellar mass, Mh/M, , and effective temperature. 
However, if we incorporate the best available spectroscopic 
determinations of stellar mass and effective temperatures to 
this analysis and expect them to be compatible with the as- 
teroseismological results, we can discriminate between the 
identification of modes and also on the other values for the 
characteristics of the models. 

In our opinion, the best fit found in the frame of case 
A mode identification, though being very good, {tp shows a 
very deep minimum) occurs however at a temperature and 
mass values too low compared to those determined by KA. 
Nevertheless, for case B we find also a very deep minimum 
at an effective temperature compatible to that of KA and a 
mass value almost identical to that given in KA. If we accept 
that this coincidence is not by chance, we have to conclude 
that this is the best fit to the structure and evolutionary 
conditions for G117B-15A. Thus, our results strongly sug- 
gest that the hydrogen envelope present in G117B-15A is a 
thick one with log (Mh/M,) ^ -3.8. 

The result regarding the actual amount of hydrogen 
present in G117B-15A should be considered as a very in- 
teresting one. In fact this is a very thick hydrogen envelope. 
Notice that this is the upper limit we considered for this 
quantity. Because of this fact, in improving our fit we may 
have considered even thicker envelopes (see Figs. 4-6). How- 
ever, notice that performing asteroseismology in the way we 
have done, it is not clear if such a procedure would produce 
stellar models plausible from a physical point of view. Notice 
that even in the hypothetical case that our artificial proce- 
dure for generating the stellar models were able to handle 
such models, there is no way to be confident that previous 
nuclear burning were effective enough to prevent the exis- 
tence of such structures in Nature. Because of this reason 
we think it would be a more safe procedure to compute the 
structure of objects of the mass of G117B-15A starting from 
the hydrogen burning main sequence throughout all the evo- 
lutionary stages previous to the WD stage. This is the only 
way we envisage in order to assess the possible existence of 
hydrogen envelopes thicker than those we have considered 
here. In any case, we should point out that previous detailed 
calculations of the evolution of objects that finally form a 
WD predict values near Mh = 1O"*M0 for 0.60 Mgobjects 
at the start of the cooling branch. However a serious study 
of this problem would carry us too far afield and is beyond 
the scope of the present paper. 

Finally, we can revisit the estimation of the parallax of 
G117B-15A as done in Bradley (1998). Because the effective 
temperature is higher (11790K vs. 11600K) and the stellar 
mass lower (0.52 Mqvs. 0.60 Mq, and thus a larger radius), 
the luminosity (log (L/Lq) =-2.44 vs. log (L/Lq) =-2.52) is 
higher than predicted by Bradley (1998) and for the same 
visual apparent magnitude mv ~ 11.50 and bolometric cor- 
rection of -0.611 mag (Bergeron, Wesemael & Beauchamp 
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1995) we find that G117B-15A should have a parallax of 
15.89 mas whereas Bradley determined it to be 16.5, which 
is still too high to be compatible (to a la level) with the 
determination of 10.5 ± 4.2 mas by van Altena, Lee & Hof- 
fleit (1994). In view of the amount of compatible evidence 
on the value of the stellar mass and effective temperature 
we expect the observed parallax of G117B-15A to be under- 
estimated. It seems difficult to adjudicate the discrepancy 
in parallax as due to an overestimation of the value of the 
stellar mass. Were this the case, it should be significantly 
lower than found in this work. 

Table 2. Summary of results for the best fitting 



Quantity 


Fitted model 


G117-B15A 




0.525 


0.53 




7.85 


7.86 ± 0.14 


Teff [K] 


11,790 


11,900 ± 140 


log(MH/M.) 


-3.83 




log(MHo/M,) 


-2.00 


10.5 ± 4.2 


Parallax [mas] 


15.89 


P{k = 2) [s] 


208.80 


215.2 


P{k = 3) 


278.85 


271.0 


P(fc = 4) 


308.57 


304.4 


P{k = 2) [10-^^s s"^] 


4.43 


2.3 ± 1.4 


P{k = 3) 


3.22 




P{k = 4) 


5.76 




References: (1) Koester & AUard (2000), (2) van Altena et al. 
(1994), (3) Kepler et al. (1982), (4) Kepler et al. (2000) 



6 CONCLUSIONS 

In this paper we have studied the structural characteristic of 
the variable DA white dwarf (WD) G117B-15A by applying 
the methods of asteroseismology. In doing so we have con- 
structed models of WD evolution considering updated and 
detailed physical ingredients. It should be remarked than we 
have included several processes responsible for the diffusion 
of elements in the WD interior. In particular, we consid- 
ered gravitational settling, chemical and thermal diffusion. 
Starting from an artificial model, we have considered several 
thickness for the outermost hydrogen layer, whereas for the 
inner helium-, carbon- and oxygen-rich layer we considered 
the profiles predicted by Salaris et al. (1997). The range of 
stellar masses we have considered is 0.50 < M^/Mq < 0.60. 
As far as we are aware, this is the first study in which evo- 
lutionary models of DA WDs considering element diffusion 
have been constructed for a set of values of hydrogen mass 
fraction. 

The evolution of each model sequences were followed 
down to effective temperature of 12500K from where on, 
we considered the evolution coupled to the oscillations. 
We considered dipolar, adiabatic g-modes with radial order 
k = 1, ■ • ■ , 4. After constructing the full set of models, we 
considered the location of the minima of a function conve- 
niently defined for the purpose of fitting. We considered two 
possible mode identifications according to Bradley (1998): 
that observed modes 215.2, 271 and 304.4 s correspond to 



k — 1, 2, 3 and fc = 2, 3, 4. We find that the employment of 
asteroseismology does not provide an univocal answer about 
the correct identification of modes present in G117B-15A. 
However, if we use them together with data deduced from 
spectroscopy it is found that the identification of the ob- 
served periods of 215.2, 271 and 304.4 s with dipolar modes 
with fc = 2, 3, 4 is clearly the best. 

The results presented in this work (summarised in 
table 2) strongly suggest that G117B-15A is a DA WD 
with a stellar mass of 0.525 Mq, a hydrogen mass 
fraction log(MH/M«) =-3.83 and effective temperature 
Tefi «11800K. Notably, the values of the effective tempera- 
ture and stellar mass are in very nice agreement with those 
predicted by spectroscopic analysis by Koester & AUard 
(2000). This represents the main result of the present pa- 
per. 

While the favoured value for the mass fraction of hydro- 
gen is the maximum we have considered, we have not tried to 
study the case of slightly (25% say) thicker hydrogen layers 
because our artificial starting evolutionary technique pre- 
vents us from being confident with the possibility that such 
structures can be actually be formed in Nature. The only 
way to answer this question is to perform full evolutionary 
calculations from the initial stages of evolution (hydrogen 
main sequence) to the WD stage. 

Thus, performing full evolutionary studies coupled to 
oscillations is, in principle, interesting in two senses. First, 
it can provide significative improvements to the fitting of 
the oscillatory modes we have found for G117B-15A in this 
paper. Second, and more importantly, it can provide us with 
a stringent test of the quality of detailed theoretical evolu- 
tionary models of WD stars by comparing the theoretical 
pattern of non-radial oscillations against accurate observa- 
tional determinations of the oscillatory pattern of WD stars. 
Work in this sense is in progress and will be published else- 
where. 
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Figure 1. The internal cliemical profile of a 0.55 Mq WD model for ^H, ^He, "He, ^^c, 16q^ £qj. ^^le case of log(MH/M,) 

=-3.862 at logT(,ff =4.096. For more details, see text. 
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Figure 2. The profile of the bottom of the hydrogen envelope for our suite of 0.55 AfQmodels at T^ff =12,500K represented by solid 
lines. For the corresponding values of Mh/M» , see Table 1. Thin lines represent the profiles corresponding to the standard treatment of 
equilibrium diffusion in the trace element approximation for the same Mh/M* values. Notice that, while for thick hydrogen envelopes the 
profiles are rather similar, there appear large differences in the case of thin envelopes. Circles denote the change of slope in the profiles 
calculated with equilibrium diffusion. For more details, see text. 
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Figure 3. Some of the main characteristics of our models for the case of a mass of 0.55 Mq. Upper panel: carbon (solid hnes) and 
helium (long dashed lines) profiles. Middle panel: The Ledoux term of the Brunt- Vaisala frequency. Lower panel: The logarithm 
of the squared Brunt- Vaisala frequency. From the centre outwards, the three peaks of the Ledoux term B are due to steep slopes in 
the carbon profile. The first two are due to the structure of the carbon - oxygen interface whereas the other is due to the helium - 
carbon interface. The inset in lower panel shows the behaviour of the Brunt- Vaisala frequency in the case of time dependent diffusion 
and according to the predictions of equilibrium diffusion for our sequence with the thickest hydrogen envelope at the hydrogen - helium 
interface (thick and thin lines, respectively). For more details, see text. 
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Figure 7. The value of ?/) as a function of the effective temperature of the star assuming case B for the largest values of log (Mh/M,) . 
Notice that if we assume the effective temperature to be in the range of the determination of Koester & AUard (2000), then, we naturally 
find that the mass of the object should be 0.525 Mq and its hydrogen mass fraction should be log (Mh/M*) =-3.83. 
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Figure 8. The periods of oscillation for radial order k = 2 (lower panel) fc = 3 (middle panel) and k = 4 (upper panel) for the case 
of models with the thickest hydrogen envelope for each considered value of the stellar mass (from bottom to top, curves correspond to 
stellar masses of 0.60, 0.55, 0.5375, 0.525, 0.50 Mq). Notice that the fitting to the observed modes (denoted with horizontal short dashed 
lines) should be considered simultaneously at the same effective temperature. 
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